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Q^. Abstract. We investigate constraints on dark energy fluctuations using type la 

^ I supernovae. If dark energy is not in the form of a cosmological constant, that is 

\—( • if the equation of state w ^ —1, we expect not only temporal, but also spatial 

C/3 , variations in the energy density. Such fluctuations would cause local variations in 

^ ' the universal expansion rate and directional dependences in the redshift-distance 

relation. We present a scheme for relating a power spectrum of dark energy fluctuations 
I to an angular covariance function of standard candle magnitude fluctuations. The 

' predictions for a phenomenological model of dark energy fluctuations are compared 

CO , to observational data in the form of the measured angular covariance of Hubble 

diagram magnitude residuals for type la supernovae in the Union2 compilation. The 
observational result is consistent with zero dark energy fluctuations. However, due 
to the limitations in statistics, current data still allow for quite general dark energy 



' fluctuations as long as they are in the linear regime 
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. y3 , 1. Introduction 

The current standard model of cosmology has the following, at least approximate, 
properties [H El E] : 

• It is spatially fiat, i.e., the total density is close to the critical density. 

• 5 % of this density is made up of baryonic matter. 

• 25 % is made up of pressureless dark matter, so far only observed through 
its gravitational effect on the universal expansion, structure formation and the 
dynamics of collapsed structures down to the scale of dwarf galaxies. 

• 70 % in a cosmological constant, observed through an apparent accelerated universal 
expansion at low redshifts. 
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However, there is still room for sizeable deviations from the concordance picture. For 
example, current data allow for the dominant energy contribution to have properties 
quite different from that of a cosmo logical constant. Although the equation of state, 
relating the pressure and the density by a linear relation, p = wp, is probably close 
to the cosmological constant value w = —1 at low redshifts, it may have very different 
properties at higher redshifts [1]. 

A dark energy component with w —1 not only allows for a time varying energy 
density, but will also exhibit spatial fluctuations. In recent years, several studies have 
been carried out concerning clustering properties of different dark energy models (see, 
e.g., Mota et al [4J and references therein). Depending on the dark energy model 
employed, clustering can occur on very different scales and at different amplitudes 
and have different couplings to the dark matter clustering. As a concrete example 
of clustering scales and amplitudes in the non-linear regime, it was shown in Mota et al 
[4j that the density contrast of a dark energy scalar field can be of order 10~^ (1 + w) in 
voids of radii 100 — 300 Mpc. The clustering can be anti-correlated with the dark matter 
clustering, which is the case for, e.g., a quintessence field with a quadratic potential in 
the linear regime [5]. Clustering of dark energy may also affect structure growth in 
the dark matter component by, e.g., introducing modifications to the spherical collapse 
model [6]. In short, there are no definite theoretical predictions on how dark energy will 
cluster in the universe. Coupled to the fact that it is notoriously difficult to measure 
temporal variations in the dark energy density, this warrants observational constraints 
on the amount of spatial clustering of dark energy in order to rule out or to further 
confirm the cosmological constant hypothesis. 

Dark energy clustering will manifest itself through several different observational 
effects. Dark energy inhomogeneities will change the form and time evolution of 
gravitational potentials to affect cosmic microwave background (CMB) anisotropies 
and the late-time integrated Sachs- Wolfe (ISW) effect, structure growth will be 
modified, and the properties of collapsed halos can change to affect cluster number 
counts [aiHl[9l[Tni[IIl[12ltt3l[Tl[T5l[l6l[I71[l8l[Tg 

In this paper we consider how dark energy clustering can affect observations of 
type la supernovae (SNe la). While the ISW effect takes into account the impact local 
inhomogeneities have on photons arriving from the last scattering surface at redshift 
z ^ 1090, observations of SNe la probe light propagation out to z ~ 1.4. When 
light from a SN la passes through a region where the dark energy density fluctuates, 
the inferred luminosity distance will change due to the local expansion perturbation 
that the fluctuation creates. This will show up as correlated fluctuations in the SN la 
peak magnitudes, where the magnitude correlation amplitudes and angular scales will 
depend on the specific dark energy clustering amplitudes and scales. A framework for 
constraining dark energy clustering models with SN la data was considered in Cooray et 
al [20j. Observational constraints on the amount of correlation in the magnitude 
fluctuations were established in Blomqvist et al [21j. Directional dependence in the 
SN data can also be used to probe models with an anisotropic dark energy equation of 
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state [221 E3] or to search for deviations from an isotropic expansion rate [2U [251 [2Sl [27] . 
An anisotropic relation between the redshift and magnitude of SNe la is also expected 
in models where an observer is located off-center in a spherically symmetric void in the 
matter distribution [28], [29] . 

Detecting correlations in SN la data is demanding because of the scatter in the peak 
magnitudes, arising both from intrinsic variation and from systematic effects that could 
induce correlated fluctuations. The observed peak magnitude is a derived quantity that 
depends sensitively on the details of how the SN la luminosities are standardised [30] . 
Also, a number of physical effects may contaminate the results. At low redshift, 
correlated fluctuations can arise from large-scale peculiar motions [SH [32]. At higher 
redshift and small angular separations, the dominating physical effect is gravitational 
lensing [33]. Intervening galactic dust may also induce correlations on small angular 
scales [M] . 

In this paper, we consider the framework presented in Cooray et al ^20j and 
continue the development of the necessary formalism for constraining dark energy spatial 
fluctuations. We also point out some of the limitations that the formalism carries. 

The paper is organised as follows: In Section [2], we describe the method to analyse 
SN la magnitude residuals for angular correlations and apply it to current data to 
establish observational constraints on the amount of correlation. In Section[3l we develop 
a theoretical formalism to relate a phenomenological power spectrum of dark energy 
spatial fluctuations to an angular covariance in the magnitude fluctuations. Our main 
results are presented in Section [H where we compare the angular covariance predicted 
from different models to the observational constraints. This paper is concluded in 
Section [5l 

2. Angular correlations in supernova data 

In this section, we investigate whether there are correlated fluctuations in the observed 
peak magnitudes of SNe la at different angular separations. We flrst give a general 
description of how we calculate the angular covariance function and then determine the 
observational constraints on the amount of correlation over the entire sky using current 
SN la data. Throughout the paper we will assume a flat universe. 

2.1. Method 

The luminosity distance to a SN la at redshift 2; in a spatially flat, homogeneous universe 
is given by 



z 




(1) 







where 




(2) 
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Here we have assumed a generic dark energy component with an equation of state (z) 
whose energy density evolves with redshift as 



l + z' 



(3) 



The apparent magnitude m of the SN la is related to the luminosity distance as 

where M is the absolute magnitude. 

We apply a methodology similar to that outlined in Blomqvist et al ^21j to 
investigate angular correlations in SN la data. For all SNe la we have redshifts, z, 
observed peak magnitudes, mobs, with uncertainties, cXm,, as well as positions in the sky. 
We form all possible SN la pairs and calculate their angular separations, 9, in the sky. 

The magnitude residual, (5m, of a SN la on the Hubble diagram is the difference 
between its observed peak magnitude and the magnitude predicted from its redshift in 
the best fit cosmology, rngt, 

5m = mobs - "ifit . (5) 
The angular (auto)covariance function as a function of angular separation is defined 

as 

cov{e) = {{5m - nsm) ■ {5m{e) - nsm)) ■ (6) 

The covariance function gives a statistical measure of the amount of correlation of 
the magnitude residuals at different angular separations. If cov{9) > 0, then SNe la 
separated by the angle 9 on average have correlated magnitude residuals, if cov {9) < 0, 
the magnitude residuals are anti-correlated. Uncorrelated magnitude residuals will give 
cov{9) = 0. Since the uncertainties 0"^ vary between SNe la, we weigh each magnitude 
residual using the weights in the calculations, fism is thus the weighted mean of 

all the magnitude residuals and the outer brackets in equation (|6]) denote a weighted 
expectation value. Writing 5m' = 5m — fism, we have explicitly that 

(^5m^5m^\ 1 ^ 

cov{9)=iY: ^2,^2 . ■ L^^^ ' (7) 

where the summation is over all SN pairs with angular separation 9. Using this 
procedure, we expect the dispersion in the covariance based on independent SN pairs 
with a given angular separation to be 

1 

' 2 



y i j m,t m,j J 

For an ideal set of SNe la with equal uncertainties Om-, the dispersion would reduce to 
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Unlon2 - 557 SNe 




Figure 1. Sky distribution in galactic coordinates of the 557 SNe la in the Union2 
data set. SNe la with z < 0.1 are marked with pluses and z > 0.1 with squares. The 
distinct feature in the lower left part consists of SNe la along the equatorial plane 
discovered by the SDSS-II survey. 

where A^p is the number of unique pairs with the given angular separation. Since each 
SN is included in many pairs, the actual dispersion will be slightly less than that given 
in equation (j8]). 

2.2. Data set 

For the data analysis, we employ the Union2 SN la data set presented in Amanullah et 
al [Ij. The compiled data set consists of 557 SNe la covering the redshift range 
z = [0.015, 1.4] and augments the Union compilation [35] by adding SNe la at low 
and intermediate z discovered by the CfA3 [36] and SDSS-II Supernova Search [37j . 
respectively, as well as six new SNe la discovered by the Hubble Space Telescope 
at high z. The peak magnitudes and uncertainties were obtained using the SALT2 
light-curve fitter. The uncertainties include both the observational and the intrinsic 
magnitude scatter. Figure [1] shows the sky distribution of the SNe la in the Union2 
data set. SNe la with z < 0.1 are marked with pluses and z > 0.1 with squares. The 
sky distribution affects on what scales and to what level we can constrain the angular 
covariance function. Whereas the low z SNe la are distributed quite evenly across the 
sky, the SNe la discovered at high z are mainly confined to small survey patches. The 
most distinct feature in Figure [1] is the dense collection of SNe la along the equatorial 
plane discovered by the SDSS-II survey. We have augmented the Union2 data by adding 
the SN positions. The data used in the analysis are presented in Table [^. 

I When used, please cite Amanullah et al [1] in addition to this paper. The Union2 data, including 
the covariance matrix with systematics, are available at |http : / / supernova . Ibl . gov/Union 



Constraining dark energy fluctuations with supernova correlations 



6 



SN 


z 


M 




RA(J2000.0) 


Dec(J2000.0) 


1993ah 


0.0285 


35.34 


0.23 


23:51:50.200 


-27:57:42.00 


1993ag 


0.0500 


36.68 


0.17 


10:03:34.500 


-35:27:45.00 


19930 


0.0529 


36.82 


0.16 


13:31:08.370 


-33:12:54.60 


1993b 


0.0701 


37.44 


0.16 


10:34:52.180 


-34:26:30.20 


1992bs 


0.0627 


37.48 


0.16 


03:29:29.970 


-37:16:18.60 



Table 1. Union2 compilation. Redshifts are given in a frame at rest with respect to 
the CMB. The distance modulus is defined as /i = to — Af. This table is available in its 
entirety in a machine- readable form at http : // www, astro . su. se/ -michaelb/SNdata^ 
A portion is shown here for guidance regarding its form and content. 

2.3. Results 

We fit a spatially flat model with dark energy with a constant equation of state and 
cold dark matter {wCDM model) to the data using the full covariance matrix with 
systematic errors. The best fit cosmological parameters are = —1.2 and fi^ = 0.34. 
The magnitude residuals are then obtained from equation ([5]). 

The angular covariance function for the 557 SNe la in the Union2 data set is 
presented in Figure O The data have been divided into ten uniform bins giving an 
angular resolution of 18°. Note that since each SN is included in several bins, the data 
points will not be independent. The error bars represent the 95% confidence limit. 
The data set is consistent with the SN magnitude residuals being uncorrelated with 
the covariance constrained to < 10^^ over the entire sky. The error bars are obtained 
using a bootstrap resampling method where we generate new data sets by keeping the 
positions in the sky fixed and randomly picking one of the SNe la for each position. The 
covariance function is then calculated for 10 000 data sets to estimate the spread in the 
covariance function. We point out that our results are not particularly sensitive to what 
cosmology we subtract when calculating the magnitude residuals; using instead the best 
fit fiat ACDM model has very small impact on the covariance function. In the next 
section, we investigate how we can use our observed limit on the angular correlation of 
SN la magnitude residuals to constrain dark energy fluctuations. 

3. Angular correlations from dark energy fluctuations 

In this section, we consider how spatial fluctuations in the dark energy density induce 
correlated fluctuations in the SN la magnitudes. We will employ a formalism that 
has previously been used in galaxy counts correlations and weak gravitational lensing 
analysis, and adapt it to relate the spatial power spectrum of dark energy fluctuations 
to an angular covariance function of magnitude fluctuations. Our goal is to estimate the 
possible effects that dark energy fluctuations can have on the SN covariance function 
which can be tested against current data. The formalism that we develop extends 
that presented in Cooray et al [20] but is not without restrictions. As such, we only 
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Figure 2. Angular covariance function for the 557 SNe la in the Union2 data set. The 
error bars represent the 95 % confidence hniit. The data have been binned uniformly 
using an angular resolution of 18°. The horizontal bars indicate the range of each bin 
and the points arc placed at the average angular separation in each bin. Note that the 
data points are not independent, since each SN contributes to several bins. 



consider it to be a toy model, serving as a first step towards a more complete formalism. 
In our setup we use a phenomenological dark energy power spectrum as input for our 
model, neglecting the observed matter power spectrum. The main limitation, however, 
is that it is currently not possible to calculate exact distances in general inhomogeneous 
cosmologies, and we will have to rely on simplified models for relating dark energy 
fiuctuations to luminosity distance fiuctuations as explained below. 

3.1. Power spectrum projection 

Suppose that there is a dark energy density fiuctuation field defined by 

^.(n..)^ "-'"-:'-/-'^' , (10) 

Px [Z) 

where the bar indicates spatially averaged quantities and fi denotes directional 
dependence. We assume that fluctuations in the dark energy density will perturb the 
expansion rate locally. Consider a perturbed Hubble function with perturbations in 
both the dark energy and matter component, 

SttG 

H"^ (n, z) = [pm {z) + {z) + 5px (n, z) + 5p^ (n, z)] . (11) 

The matter perturbation is introduced to balance the dark energy perturbation such that 
we recover H{z — Q) = Hq, i.e., the Hubble constant does not have spatial variations. 
Alternatively, one could choose to normalise the local expansion such that the Hubble 
parameter is spatially homogeneous at an earlier time slicing, or to impose a condition 
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of homogeneous Big Bang time for each local patch. Using this setup we can write (see 



Appendix A for more details) 

(n, z) = Hi [n^ (1 + zf + {z) + 5. (fi, z) {z) a {z)] , (12) 

where 

a{z) = l- -^^T^^ (1 + z)-""-'^^'-'^ (13) 

ensures that we have a spatially homogeneous expansion at z = 0. The parameters 
Sm and Sx determine the degree of gravitational growth with redshift of the matter 
and dark energy density contrast, respectively. We point out that we have assumed a 
flat universe, and as a consequence the perturbations of dark energy and matter are 
necessarily anti-correlated. A less restricted setup of the perturbations would also allow 
for spatially varying curvature, but unfortunately it is currently not possible to calculate 
exact distances in such a general model. 

The presence of dark energy inhomogeneities along the line of sight to a SN la 
will cause fluctuations in the observed peak magnitude, or equivalently, a fractional 
fluctuation in the luminosity distance, related by 

5m = . (14) 

In 10 dL 

The averaged fractional luminosity distance fluctuation integrated over all sources in 
the direction fi, out to the most distant SN la at comoving distance r^, is given by a 
projection of the dark energy fluctuation field down on the celestial sphere, 

-l-(n) 



average ^ 



drW{r)6^[n,z{r)] . (15) 
The exact nature of the projection is determined by the weight function W (r), derived 



m 



Appendix A 



W (r) = -^^^^/ (r)] a [z (r)] J dr'^ . (16) 

r 

The weight function basically tells us how much a dark energy fluctuation at a comoving 
distance r will affect the luminosity distance of SNe at larger r. A key assumption in 
the derivation of the weight function is that the density fluctuations are small, i.e., dark 
energy clustering only occurs in the linear regime. Here n (r) is the radial distribution 
of SNe la, normalised so that 



drn (r) = 1 . (17) 



To calculate the induced correlations in SN la magnitude fluctuations it is enough 
to know the statistical properties of the fluctuation field. If the dark energy density 
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contrast is a Gaussian random field, these are encoded in the spatial power spectrum 
Px (k), formally defined as 

(5. (k) 6: (k')) = {2nf Si (k - k') Px (k) . (18) 

Here k is the comoving wavenumber related to the comoving wavelength Xhj k = 2'k/X. 
The spatial power spectrum of dark energy fluctuations can be related to an angular 
power spectrum of fractional luminosity distance fluctuations Ci using a projection in 
the flat-sky approximation involving the same weight function, 

Q= dr'^P^[k=lz{r)^ . (19) 



The angular power spectrum tells us how much fluctuations there are on angular scales 
6 = 180°//. Equation f ll9p is a projection down on a small patch of the sky that 
effectively can be considered flat, and is only valid for small angular separations. This 
form of the flat-sky angular power spectrum requires that the weight function W (r) 
is fairly constant over the regions of interest and that the correlation lengths are small 
compared to the integration range. This relation has been used extensively in relation to 
galaxy counts [38] and weak lensing observations [391 SO]- We note that a more general 
projection down on the entire celestial sphere, valid for all angular separations, is indeed 
possible to perform, but is more cumbersome to work with. Since we expect the induced 
correlations to be largest at small angular separations, equation f fT9|) is sufficient for this 
study. 

From the angular power spectrum of fractional luminosity distance fluctuations, we 
calculate the angular covariance function of magnitude fluctuations as 

^""'"'"(hno) ItJ^m'')' (2°) 



where Jo is the zeroth order Bessel function of the first kind, and the prefactor comes 
from equation f|T^ . 

In summary, we use equations (fT6l) . ( fT9l) and ( !20l) to relate the spatial power 
spectrum of dark energy fluctuations to the angular covariance function of magnitude 
fluctuations. 

3.2. Phenomenological power spectrum of dark energy fluctuations 

In order to apply our formalism to constrain dark energy clustering, we need to specify a 
power spectrum of dark energy fluctuations that we can use to compare the theoretical 
covariance from the projected power spectrum with the observational covariance from 
SN la data. As mentioned in the introduction, dark energy clustering is very model 
dependent. Thus, from a theoretical perspective, it is not possible to predict a single 
specific power spectrum. We will therefore use a phenomenological power spectrum 
Px(A;, z) which covers the amplitude distribution of different wavenumbers, with a cut-off 
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Figure 3. Phenomenological dark energy spatial power spectrum today for = 
[1.5, 1,0.5]. The left panel shows models with cut-off wavenumber fee = 0.1 h Mpc^^. 
The amplitude values are (in the above order) = [10~^^,2 x 10^^^, 10^^°]. The 
right panel instead shows models with fee = 0.05 h Mpc~^. Here the amplitudes are 
(5^jj = [10^^^, 3 X 10~^^, 10~^]. The dot-dashed line indicates the linear matter power 
spectrum from Peacock 41, . 



towards small fluctuation scales and a redsliift dependent growth factor. The specific 
parameterisation was previously considered in Cooray et al [20] and is basically a 
modified inflation power spectrum. It is given by 

(fc) ^ ^^3^SM. = ( e-'/'^ (1 + z)-- . (21) 



27r2 \Ho^ 

Here we have also introduced the dimensionless quantity A'^{k), which gives the 
contribution to the variance of the fluctuation field per logarithmic interval. The power 
spectrum has four free parameters which can be constrained by data: S^^i gives the 
fluctuation amplitude on horizon scales; gives the power law distribution of the 
fluctuations; kc is a cut-off wavenumber above which fluctuations are suppressed; and 
Sx is the growth parameter from equation f|T3l) . 

In the analysis, we fix the growth parameters to Sm = Sx = 2 corresponding to 
a growth directly proportional to the scale factor (i.e. the density contrast decreases 
with redshift) as expected for linear matter fluctuations in a matter dominated universe. 
Since we only integrate out to z = 1.4 in the power spectrum projection, the final results 
are not sensitive to this choice. 

Figure [3] shows the power spectrum for different values of the parameters and 
kc- The values considered ensure that the spectra tend to zero on large scales, where the 
projection formalism breaks down, and that fluctuations are suppressed on small scales 
where clustering is not expected. The curves have been normalised so that the variance 
of the fluctuation field is (5^) = 0.1, since we demand that the density fluctuations are 
small. This means that the amplitude 5xh been adjusted for each individual model 
according to this normalization condition (see the figure caption for the different values 
used). Also shown for comparison is the linear matter power spectrum calculated using 
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Figure 4. Weight function used in the power spectrum projection. It is based on the 
radial distribution of the 392 SNe at z > 0.1 in the Union2 data set. 



a fitting formula in Peacock [31] . Tlie dark energy power spectra all lie below tlie matter 
power spectrum. Tfiis need not necessarily be expected from a dark energy clustering 
model, but is instead only a consequence of imposing that the fluctuations are small in 
our formalism. 

Given the power spectrum of dark energy fluctuations, we can calculate the angular 
power spectrum of fractional luminosity distance fluctuations using equation (fT9l) . Note 
that the angular power spectrum depends on the redshift distribution of SNe la through 
the weight function. The radial SN distribution n(r) is taken directly from the Union2 
data set. The distribution is thus a sum of delta functions that pick out the specific 
distances to the different SNe la such that the weight function W (r) becomes 



where A'^sn is the number of SNe la in the data set and the summation is over all 
distances larger than r out to the most distant SN at rg. As the unperturbed background 
cosmology, we use the best fit fiat wCDM model in Section 12.31 

We choose to make a redshift cut in the data set, excluding all SNe la at 
z < 0.1, for a couple of reasons. Firstly, the non-uniform radial distribution of 
SNe la, which is particularly populated at low z, gives a spike in the weight function 
at r = 50 — 100 h^^ Mpc. The redshift cut removes this feature, giving a smoothly 
varying weight function as necessary in the power spectrum projection. Secondly, low 
z SNe are affected by their peculiar motions which can induce correlated fluctuations 
in the observed peak magnitudes. Removing all SNe at z < 0.1 is a conservative choice 
to obviate this effect. The redshift cut brings down the total number of SNe la in the 
data set from 557 to 392 (see also Figure [1] for which SNe that are removed). The 
observational limits used in the next section will essentially be unaffected by this cut. 




(22) 
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Figure 5. Angular power spectra of fractional luminosity distance fluctuations for 
different values of the parameters rix and kc ■ The curves are normalised as in Figure |31 



Figure H] shows the weight function in equation fl22|) based on the radial distribution of 
the 392 SNe at 2; > 0.1 in the Union2 data set. 

Figure [5] shows the angular power spectrum of fractional luminosity distance 
fluctuations. The amplitude of the curves are normalised according to the values in 
Figure [31 It would have been interesting to compare this result to the angular power 
spectrum predicted in a model with perturbations in either dark energy or matter. 
However, the current formalism considers perturbations in both energy components and, 
as such, cannot address, e.g., luminosity distance fluctuations from matter perturbations 
in a ACDM universe. 

4. Constraints on dark energy spatial fluctuations 

Using equation ( l20l) . we calculate the angular covariance function of magnitude 
fluctuations based on the 392 SNe la at 2; > 0.1 in the Union2 data set. Figure [6] 
shows the results for different values of the parameters and kc together with the 
observational limits established from the SN la data. The amplitudes of the curves 
are normalised according to the values in Figure |3l Note that we have zoomed in on 
small angular separations, since the angular covariance function increases as the angular 
separation decreases, and that the observational limits are based on a single data bin 
out to 6' = 1.5° in which there are 2086 SN pairs. We point out that the observational 
limits are essentially unaffected by the low-redshift cut made in the data set, since the 
small angle bin is vastly dominated by SN pairs at higher z. Also plotted is the expected 
angular covariance from gravitational lensing based on the linear matter power spectrum 
calculated following Cooray et al [33]. The lensing curve lies below the dark energy 
clustering curves, but we expect it to increase at the smallest angular separations if the 
non-linear matter clustering contribution would be included. In fact, there are tentative 
detections of the gravitational lensing of SNe la from non-linear matter clustering from 
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Figure 6. Angular covariance function of magnitude fluctuations for different values 
of the parameters and fee of the dark energy power spectrum. The curves 
are normalised as in Figure [3] The confidence limits obtained from the data are 
95% (orange) and 99.7% (yellow). Included is also the covariance expected from 
gravitational lensing (dot-dashed). 



the correlation of SN la magnitude residuals and the positions and masses of foreground 
galaxies p|li3]. 

The covariance induced by the dark energy fluctuations is within the observational 
limits and it is thus not possible to rule out any phenomenological model with linear 
dark energy fluctuations using the current SN la data. However, given more SN data, 
this method can yield important constraints also for linear fluctuations. Based on 
Figure [6], we see that we will be able to constrain linear dark energy fluctuations if 
the observational limits are reduced by an order of magnitude. Since the confidence 
limits scale roughly as ~ 1/ ^/N^, we expect this to be possible if the number of SN 
pairs in the bin increases by a factor of 100. This corresponds to having roughly ten 
times more SNe la with angular separations 6 < 1.5° compared to what is currently 
available. 

5. Conclusions and discussion 

Much of the current efforts in dark energy research are focussed on determining whether 
dark energy can be described by a cosmological constant or if it is a dynamical quantity. 
If we are to unveil the nature of dark energy, it is important to consider the possibility 
of both temporal and spatial variations. 

In this paper, we have investigated the effect of dark energy fluctuations on the 
angular covariance function of SN la magnitude residuals. In the first part, we performed 
a quantitative analysis of the covariance of magnitude residuals using the Union2 data 
set. We showed that the data are consistent with being uncorrelated and that the 
covariance is constrained to < 10^^ over all angular separations. In the second part, 
we constructed a formalism which relates a phenomenological dark energy fluctuation 
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power spectrum to the angular covariance of magnitude fluctuations. We found that 
the covariance induced by linear dark energy fluctuations cannot be ruled out by the 
current data. However, as the number of well observed SNe la increases, it will be 
possible to discriminate between different phenomenological models. Ten times more 
SNe with small angle separations will suffice to rule out certain models also with small 
dark energy fluctuations. Larger statistics would also allow for an analysis in terms of 
correlations in the Hubble diagram residuals as a function of spatial distances between 
SNe rather than only angular separations. 

At this point, it is important to again point out the fact that the results are not 
exact in the sense that it is not (yet) possible to calculate cosmological distances in 
general inhomogeneous cosmologies. This is usually referred to as the backreaction 
problem, i.e., how local inhomogeneities backreact on the global expansion |44]. In 
our setup, inhomogeneities produce perturbations in the Hubble function with a fixed 
background expansion. In reality, these inhomogeneities might backreact, due to the 
non-linear structure of the field equations, to produce a different global expansion rate. 
In order to obtain what we believe to be at least order of magnitude correct estimates of 
SN magnitude angular covariance, we have assumed that local fluctuations in the dark 
energy density will cause local fluctuations in the Hubble expansion that in turn will 
induce changes in the luminosity distance to a given redshift. In doing this, we have 
normalised the local expansion such that H [z = 0) = Hq, i.e., the Hubble constant does 
not have spatial variations. This is accomplished by introducing perturbations in the 
matter component by hand in the Hubble function such that the dark energy clustering 
is anti-correlated with the matter clustering. In a more refined formalism, one would 
like to include the matter and dark energy power spectrum separately. This would, 
however, require the inclusion of the simultaneous effect from local perturbations in the 
matter and dark energy component on the Hubble function. We also note that in our 
model, even for vanishing dark energy fluctuations, we would still expect an induced 
covariance from the local perturbations in the expansion history produced by the dark 
matter fluctuations. Thus, in a more general context, the method presented in this 
paper is useful for investigating Hubble flow perturbations. 

Moreover, our approach is not valid in the non-linear clustering regime, and as such 
we cannot test phenomenological models with large density fluctuations. However, since 
we expect the problem of backreaction to be more severe for non-linear fluctuations, we 
have considered it superfluous to include the extra complications of allowing for large 
fluctuations in our simple model. The formalism also breaks down at close to horizon 
scales, where dark energy fluctuations could occur. In view of these limitations, we 
welcome all efforts to solve the problem of the expansion history and light propagation 
in an inhomogeneous universe. An improved formalism will definitely be warranted in 
the event that correlations are detected in future data. 
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Appendix A. Derivation of the weight function 

We want to construct a weight function W (r) that connects the dark energy fluctuations 
(5x with the fractional fluctuations in the luminosity distance Sdj^/dj^. The goal is to arrive 
at a form 



r. 

It*"' 



j drW (r)S^[&,z(r)\ . (A.l) 



average 

^ 

We start with the Hubble function in some given direction fi (we suppress the directional 
dependence in the following equations) and add a perturbation in both the matter and 
dark energy component, 

{z) = ^ [pra (z) + Px (z) + 5p, (z) + 6pm (z)] . (A. 2) 

We normalise our model such that the Hubble constant is spatially homogeneous, 
H {z = 0) = Hq, imposing 5px,o = —Spm,o- Subscript denotes values at 2; = 0. Suppose 
that 

^ = ^ (1 + z)-/^ , (A.3) 

The parameter determines the redshift dependence of the dark energy density 
contrast. For the matter perturbation we have 



l + zy'''^'" . (A.4) 

Pm [Z) Pm,0 

Using equation (1A.3I) in equation ( ]A.4I) then gives 

6p^ (z) = -6p^ (z) (1 + . (A.5) 
The perturbed Hubble function can then be written as 

(z) = ^ [p^ (z) + px (z) + 6p^ (z) a (z)] , (A.6) 



where we have defined 



'1 + Z)^ ,^ , ,_.5^/2+s^/2 



« (z) = 1 - (1 + z)-'-^'^'-i' . (A.7) 

Dividing equation f lA.6l) by the critical density pcrit.o = 'iHl/^r^G gives 

{z) = Hi [n^ (1 + zf + fix/ {z) + 5fix {z) a {z)] , (A.8) 
where we have defined 

dil^ [z) = = _ = dx [z) iZx/ [z] . (A.9) 

Pcrit,0 P^[Z) Pcrit.O 
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We can write the perturbed Hubble function as the unperturbed Hubble function plus 
a small perturbation, 

H{z) = H (z) + 6H (z) = HoE (z) + SH (z) . (A.IO) 

Since the perturbation is assumed to be small, we can Taylor expand equation ( lA.SP to 
first order. 



H{z) = H,E {z) 

From this we identify 

5H {z) iz) 



1 + 



2^2 



-a {z) 



a {z) . 



H {z) 2^2 {z) 
For a perturbation in the direction fi, we have 



di, + 5di^ (fi) = c{l + z) 



dz' 



H {z') + 5H (fi, z') 



(A.ll) 



(A.12) 



(A.13) 



Notice that we have explicitly written the directional dependence on the perturbation 
in both the luminosity distance and the Hubble function, since the expansion history 
can be different in different directions. Doing a Taylor expansion of the integrand, we 
can identify 



Ml (n) = -cil + z) 



Shifting to the r coordinate through dr 
equation flA.121) and equation flA.9|) . we obtain 



dz' 5H{h,z') 
H {z') H {z') ■ 

[c/H{z)'j dz and using c^l 



6di^ (fi,r) 
dL (r) 



dr' 



2^2 



^f{r')a{r')5An,r') 



(A.14) 

[l + z)r, 

(A.15) 



Note that, since a{z) < for z > 0, a local dark energy overdensity will increase 
the luminosity distance in the direction of the fluctuation. We can obtain the total 
averaged fractional fluctuation in the luminosity distance integrated over all sources in 
the direction fi using the normalised source distribution function n(r). 



6di 
di, 



n 



average 



-Jdr 





n (r) 



dr' 



n 



2^2 (r 



^f{r')a{r')6An,r') .(A.16) 



Since the distribution function n (r) will be a sum of delta functions which represents the 
distribution of the SNe la, we can switch the integration limits. What equation ( ]A.16P 
essentially does, is to pick out each SN and look at the perturbations 6x (fi, r) in front of 
it to see how they affect the luminosity distance. This will be equivalent to looking at the 
perturbation 6^ (fi, r) at each point r and see how each SN behind it in the line of sight 
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will be affected in terms of the luminosity distance. We can thus write equation (IA.16P 

as 



— dr 



average ^ 



5,{n,r) .(A. 17) 



Comparing this expression with equation ( lA.ip . we can identify the weight function 
W (r) as 

W (r) = -j^^^f [z (r)] a [z (r)] j dr' ^ . (A.18) 



r 
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